Photoassociation of a Bose-Einstein Condensate near a Feshbach Resonance 



00 

o 
o 

(N 



M. Junker, D. Dries, C. Wclford, J. Hitchcock, Y. P. Chen and R. G. Hulet 

Department of Physics and Astronomy and Rice Quantum Institute, Rice University, Houston, TX 77251, USA 

(Dated: July 9, 2008) 

We measure the effect of a magnetic Feshbach resonance (FR) on the rate and light-induced 
frequency shift of a photoassociation resonance in ultracold 7 Li. The photoassociation-induced loss 
rate coefficient, K p , depends strongly on magnetic field, varying by more than a factor of 10 4 
for fields near the FR. At sufficiently high laser intensities, K p for a thermal gas decreases with 
increasing intensity, while saturation is observed for the first time in a Bose-Einstein condensate. 
The frequency shift is also strongly field-dependent and exhibits an anomalous blue-shift for fields 
just below the FR. 

PACS numbers: 03.75.Nt, 34.50.Rk, 33.80.Gj, 34.50.-s 
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Increasing interest in producing ultracold molecules 
by associating ultracold atoms has motivated an im- 
proved fundamental understanding of association pro- 
cesses. Photoassociation (PA) and magnetically-tuned 
Feshbach resonances (FR) are two ways in which gases 
of ultracold atoms may be connected to the molecular 
bound states of their underlying two-body interaction 
potentials. A Feshbach resonance is realized by tuning 
a molecular bound state near the scattering threshold 
[H- The bound state perturbs the scatteringstate which 
can, in turn, strongly affect the rate of PA [2j ■ The rate 
of association of a Bose-Einstein condensate (BEC) is 
of particular interest, as a universal regime is predicted 
at extremely high rates of association [J, 2, where 
the rate is limited by a universal transient response of the 
atom pairs 0] , independent of the underlying microscopic 
association process. 

Several groups have measured absolute PA rate coef- 
ficients HH, Ell U, 12 [. Of these, however, only two 
were performed with a quantum gas. In ref. [8[ , the rate 
of PA of a sodium BEC was found to increase linearly 
with laser intensity, without any indication of saturation. 
We previously investigated PA of 7 Li at the critical tem- 
perature for condensation T c (negligibly small conden- 
sate fraction), and observed saturation at a value con- 
sistent with the maximum two-body collision rate im- 
posed by quantum mechanical unitarity 0]. Other ex- 
periments u sing thermal gases have made similar obser- 
vations [HI, Il3l ] The predicted universal rate limit re- 
mains elusive. Previous experiments have also measured 
an intensity-dependent spectral red shift of the PA reso- 
nance [H, 0, [H| , which has been ascribed to light-induced 



coupling to the continuum [15l . Il6l . 1 1 71 ] - In this paper, 
measurements of the PA rate coefficient, K p , and the 
spectral shift are presented in the vicinity of a magnetic 
FR of 7 Li for both a thermal gas and a BEC. For fields 
near the FR, the measured values of K p vary by more 
than a factor of 10 4 . By exploiting this huge PA rate 
enhancement, the first evidence for saturation of the rate 
of PA-induced loss in a condensate is obtained. Further- 
more, the FR strongly affects both the magnitude and 



sign of the spectral shift. 

Atomic Li in the \F=1, mp=\) state is prepared in a 
sing le-beam optical dipole trap, as previously described 
[181 ] . A uniform magnetic field can be set between 480- 
900 G to access both sides of the FR located near 736 G 
[lit . Evaporative cooling is achieved by reducing the op- 
tical trap beam intensity. Either Bose condensates with 
no discernable thermal fraction, or thermal gases with 
T > T c , can be produced depending on the final optical 
trap depth. 

The PA laser propagates collinearly with the optical 
trap beam, and is tuned to resonance (peak loss) with 
a vibrational level v of the 1 3 S+ excited molecular state 
There are no unresolved rotational or hyperfine lev- 
els. For most of the measurements we tune to the v = 83 
level, located 60 GHz below the 2 2 S l 1 / 2 + 2 2 P 1 / 2 dissoci- 
ation limit j20|. The v ~ 83 level is a good compromise 
between strong free-bound overlap and relatively weak 
off-resonance scattering from the atomic transition. Ex- 
cited molecules created by the PA laser decay into pairs 
of energetic atoms that predominately escape the trap. 
The duration of the PA laser pulse r is adjusted at each 
intensity I and field to maintain the fractional loss of 
atoms to approximately 30%. The rise and fall times 
of the PA pulse are less than 0.3 /us. The number of re- 
maining atoms following the laser pulse is measured using 
on-resonant absorption imaging. 

The on-resonance loss-rate coefficient, K pi is defined 
by the time evolution of the density distribution: h{t, r) 
= —K p n 2 (t,r), where any possible time-dependence of 
K p during the PA pulse is averaged over. For the BEC 
data (Fig. 4), r is much less than either trap period, and 
collisional redistribution during the PA pulse can be ne- 
glected. In this case, we solve for the time-dependent 
density distribution analytically and extract K p by spa- 
tially integrating the density and matching the observed 
and calculated loss after a time r [1]. For much of the 
thermal data (Figs. 1 and 3), however, r is comparable to 
the radial trap period and redistribution is not negligible. 
Without accounting for the details of the time evolution, 
we take K p to be the average of two estimates: the first 
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FIG. 1: 

1.65 W/cm 2 , while r is adjusted between 0.07 and 270 ms to 
keep the fractional loss at 30 ± 6% (r oc 1/K P ). The fitted 
temperatures range between 9-18 /xK. The spread in initial 
peak density, 0.5 — 10 x 10 12 cm -3 , reflects a significant de- 
crease in density near the FR, which we ascribe to enhanced 
three-body loss. The axial and radial trapping frequencies are 
20 Hz and 3 kHz, respectively. A systematic uncertainty in 
K p of 25% reflects the uncertainty in the evolution of the den- 
sity (see text), as well as uncertainties in probe detuning and 
intensity. The statistical uncertainty is 20%. The measure- 
ment of the highest values of K v are saturated (Fig. 3) and 
are therefore lower limits. The inset shows K v (same scale as 
main figure) for three different excited state vibrational levels. 



assumes no redistribution, as for the BEC case, while the 
second assumes that the distribution remains Boltzmann 
throughout. Due to the small fractional loss, we find 
the difference between the two estimates to also be small 
(<25%), and account for it in the stated uncertainties. 

Figure 1 shows the field dependence of K p on both 
sides of the FR for a thermal gas. The data show a 
pronounced dip in K p at 710 G, and a maximum, which 
is more than 10 4 times larger, located between 725-750 
G. Measurements could not be made closer to the FR due 
to high inelastic collisional loss rates, presumably from a 
three-body process. A condensate could not be probed 
on the high-field side of the FR due to its instability for 
a < 0, but at a given field on the low-field side K p for a 
condensate was found to be approximately half that of a 
thermal gas, consistent with quantum statistics. 

The variation in K p can be understood in terms of 
the Franck-Condon principle and the dependence of the 
scattering state wavefunction f(R) with magnetic field 
(R is the intcrnuclear separation). PA predominately 
occurs at the Condon radius, Rc, defined as the outer 
classical turning point of the excited vibrational energy 
state. Since K p oc \f(Rc)\ 2 , a node in f(R) at Rc rc- 



ously observed 21 1, although with much lower contrast. 
For sufficiently low T that the collisional wave vector k 
satisfies k(R - o) < 1, the asymptotic scattering state 
f(R) ~ fc : / 2 (i? — a), and a node in f(R) appears at 
R ~ a (22j . The inset to Fig. 1 shows that the location 
of the K p minima shift to higher fields (larger a) with 
increasing v, consistent with the fact that Rc increases 
with v. The peak in K p near 736 G, on the other hand, is 
a manifestation of the enhancement of f(Rc) when \a\ is 
large due to the FR. The ratio of the asymptotic values 
of K p on the high and low field sides of the FR is ~3.5, 
which can be understood from the ratio of |/(i?)| 2 eval- 
uated using the values of a at the two asymptotic fields. 
Using a coupled-channel calculatio n 12311 with previously 
determined potentials for lithium |24j |. we find the scat- 
tering lengths at 600 G and 900 G to be a; ~ 8a a and 
ah — — 43 a , respectively, and Rc — 103 a for v = 83, 
where a a is the Bohr radius. With these values, the ex- 
pected ratio is (Rc—ah) 2 /(Rc~a-i) 2 — 2.4, in reasonable 
agreement with measurement. 

The effect of the FR on f(R) is also evident in the light- 
induced spectral shifts. We measured spectral shifts for 
both a thermal gas and a condensate for fields around 
the FR. The location of a PA resonance at a particular 
field and intensity / is determined by fitting the reso- 
nance lineshape to a Lorentzian. By taking resonance 
curves for several values of /, we observe that the spec- 
tral shift is linear in /, and extract the slope, E, plotted 
in Fig. 2. Far from the FR, the shift is red (E < 0), in 
agreement with previous measurements of single-photon 
PA @, i, [3]. As the FR is approached from low field, 
however, E changes sign, becoming large and positive 
just below the FR and large and negative just above it 
in a dispersive-like manner. The field where E vanishes 
(710 G), coincides with the location of the zero in K p 
shown in Fig. 1. 

The spectral shift is a consequence of light-induced 
mixing of nearby molecular states into f(R), and has 
been shown to scale linearly with f(Rc) [16| . In previ- 
ous single-photon experiments the shift was dominated 
by coupling to the continuum, resulting in a red shift. In 
the case of PA near a FR, however, the shift is strongly 
affected by the underlying closed-channel molecular state 
responsible for the FR. For fields below the FR, the con- 
tribution to the shift from the bound state is positive, 
while above resonance it is negative. 

The data of Fig. 1 demonstrate that K p can be extraor- 
dinarily large near the FR, making it an ideal system for 
exploring saturation. Figure 3 shows measured values 
of Kp vs. / for a thermal gas at a field just above the 
FR. This data can be analyzed in terms of a Fano model 
of a bound state coupled to a continuum, where the on- 
resonance intensity-dependent loss rate coefficient can be 
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FIG. 2: Slope of the intensity-dependent spectral light shift 
E. The error bars are the standard error of the fitted slope at 
each field. For most of the data, this statistical uncertainty 
is smaller than the size of the plotted points. The systematic 
uncertainty in E is 10% due to uncertainty in I. Condensates 
and thermal gases were found to exhibit the same E. 



written as [l6| 

K p (I)=4K ma I M Isat (1) 

(-1 + J- sat) 

where K max is the maximum loss rate coefficient and 
I sat is the corresponding saturation intensity that ac- 
counts for broadening of the PA resonance. Because of 
the low temperatures in our experiment we neglect any 
energy dependence of K p . The solid line in Fig. 3 is a 
fit of the data to Eq. Q] The data clearly show the pre- 
dicted rollover in K p that occurs when / > I sat . We 
also observe significant power broadening that accompa- 
nies saturation. By fitting the measured lincshape to a 
power-broadened Lorentzian at 35 W/cm 2 , for example, 
we obtain I sat = 4.3 W/cm 2 , in good agreement with the 
fit to Eq.Q] which yields 4.1 W/cm 2 . 

The theory of rcf. [HI provides an estimate of 
the unitarity-limited maximum rate coefficient K u = 
2/cbTA 3 /h, where h is Planck's constant and A is the 
thermal de Broglie wavelength using the reduced mass, 
and the factor of 2 accounts for losing both atoms per 
PA event. For T = 16 fiK, K u = 8.5 x 1CT 9 cm 3 /s, in 
excellent agreement with the observed maximum. 

Figure 4 shows the measured intensity dependence of 
K p for a BEC near the FR. By achieving extremely large 
FR-enhanced loss rates, saturation was observed in a con- 
densate for the first time. The maximum K p of 1.4 x 
1CP 7 cm 3 /s is nearly a factor of 10 larger than for any 
previous PA rate measurement 

A complete theoretical understanding of PA of a BEC 
in the high-intensity regime is still evolving, but a framc- 



FIG. 3: K p vs. / in a thermal gas at 755 G. T ~ 16 /iK and 
the peak density is 7.5 x 10 11 cm -3 , r is between 210-250 /is 
for all of the data except at the lowest intensity, for which r 
= 2 ms. The error bars are the statistical uncertainties from 
the fitted T. The systematic uncertainties are 30% for K p 
and 10% for /. The solid line is a weighted fit to the data 
using Eq. 1, with the fitted parameters K max = 9.5 x 10~ 9 
cm 3 /s and I sa t = 4.1 W/cm 2 . A linear fit to the slope for 
small / gives 8.2 ± 1.1 x 10 -9 cm 3 s _1 /( w cm -2 ), where the 
uncertainty refers only to statistical uncertainty. 



work for comparison of the data may be established by 
defining K p in terms of a characteristic length, L, as K p = 
(fi/m)L. The smallest relevant length scale is the average 
interatomic separation, evaluated at the peak density n . 

— 1/3 

Taking L = n , we obtain the "rogue photodissocia- 
tion" limit, K p d, a universal regime where saturation in 
the rate of condensate loss is predicted due to the disso- 
ciation of excited molecules into the hot pair continuum 
Q. The NIST BEC experiment achieved K p ~ K pd , 
but with no indication of saturation For the data of 
Fig. 4, n a = 1.6 x 10 12 cm" 3 , giving K pd ~ 8 x 1CT 9 
cm 3 /s. Our measured K max is nearly 20 times greater 
than K p d- More recent calculations [a, [jj, however, have 
shown that while dissociation does impose a rate limit 
on condensate loss, it is not as stringent as K p d- While 
the highest-/ data of Fig. 4 are well into the "dissocia- 
tion regime" 0, 0| , our measured K max is, nonetheless, 
nearly 7 times greater then predicted from the equations 
given in ref. Q. We note that in the experiment atoms 
removed from the condensate but left with energies below 
the trap depth (500 nK) may be indistinguishable from 
condensate atoms by our detection method. This effect 
cannot explain the discrepancy with theory, however, as 
it would result in an apparent suppression of K p at the 
onset of photodissociation. 

The observed saturation could be explained by a higher 
than expected limit imposed by dissociation, perhaps due 
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FIG. 4: K p for a BEC at 732 G (a ~ 1000 a ). The 1/e 2 
intensity radius of the PA beam of 560 fj,m is much larger 
than the radial Thomas-Fermi radius of the BEC. The pulse 
duration r oc 1/K P , is 3 fj,s for the highest / and 50 /zs for 
the lowest. The peak density is 1.6 x 10 12 cm~ 3 . The axial 
and radial trapping frequencies are 4.5 Hz and 200 Hz, re- 
spectively. The solid line is a fit to Eq. 1 with fit parameters 
K max = 1.4 x 10 -7 cm 3 /s and I sat = 26 W/cm 2 . A linear 
fit to the low intensity data gives a slope of 1.6 x 10 -8 cm 3 
s _1 /(W cm -2 ), as indicated by the dashed line. The error 
bars are the statistical uncertainties, due mainly to the mea- 
sured Thomas-Fermi radius of the initial density distribution. 
A systematic uncertainty of 12% reflects the uncertainty in 
the image probe detuning and intensity. / could not be in- 
creased further without producing loss from dipole forces. 



to cross-coupling between the PA and Feshbach reso- 
nances 25 1 . An alternative explanation is provided by 
quantum mechanical unitarity. If we take L ~ 2Rtf, 
where Rtf — 10 /im is the radial Thomas-Fermi radius, 
then K u ~ 1.8 X 10~ 7 cm 3 /s, in good agreement with 
the measured value of 1.4 x 1(K 7 cm 3 /s. A preprint has 
recently appeared that directly models our experiment 
with good quantitative results [25j |. but it docs not de- 
termine whether the rate limit is imposed by dissociation 
or unitarity. 

We have presented well-characterized measurements of 
PA near a FR, which we hope will serve to constrain 
and guide theory, especially in the previously unexplored 
strong coupling regime. Condensate loss rates as high as 
the unitarity limit have been demonstrated by combin- 
ing PA and Feshbach resonances. Future studies should 
explore the dynamical behavior of loss since K p is pre- 
dicted to be time-dependent in the dissociation regime 
B HSU- Finally, Feshbach enhanced PA may prove 
useful in realizing the long-sought goal of coherent os- 
cillation between atomic and molecular condensates 26 1 
or for efficient production of ground-state molecules 27 1 



using a two-photon process [28j. 
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